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8.1 Introduction 

8.1.1 Tunneling in integrable systems 



Since the early days of quantum mechanics, tunneling has been recognized as one of the hall- 
marks of the wave character of microscopic physics. The possibility of a quantum particle to 
penetrate an energetic barrier represents certainly one of the most spectacular implications 
of quantum theory and has lead to various applications in nuclear, atomic and molecular 
physics as well as in mesoscopic science. Typical scenarios in which tunneling manifests are 
the escape of a quantum particle from a quasi-bounded region, the transition between two 
or more symmetry-related, but classically disconnected wells (which we shall focus on in the 
following), as well as scattering or transport through potential barriers. The spectrum of 
scenarios becomes even richer when the concept of tunneling is generalized to any kind of 
classically forbidden transitions in phase space, i.e. to transitions that are not necessarily 
inhibited by static potential barriers but by some other constraints of the underlying clas- 
sical dynamics (such as integrals of motion). Such "dynamical tunneling" processes arise 
frequently in molecular systems P and were realized with cold atoms propagating in period- 
ically modulated optical lattices ^ IH H]. Moreover, the electromagnetic analog of dynamical 
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tunneling was also obtained with microwaves in billiards [5]. 

Despite its genuinely quantal nature, tunneling is strongly influenced by the structure of 
the underlying classical phase space (see Ref. [B] for a review). This is best illustrated within 
the textbook example of a one-dimensional symmetric double-well potential. In this simple 
case, the eigenvalue problem can be straightforwardly solved with the standard Jeffreys- 
Wentzel-Kramers-Brillouin (JWKB) ansatz [7]. The eigenstates of this system are, below 
the barrier height, obtained by the symmetric and antisymmetric linear combination of 
the local "quasi-modes" (i.e., of the wave functions that are semiclassically constructed on 
the quantized orbits within each well, without taking into account the classically forbidden 
coupling between the wells), and the splitting of their energies is given by an expression of 
the form 



Here E is the mean energy of the doublet, V{x) represents the double well potential, m is the 
mass of the particle, Q denotes the oscillation frequency within each well, and the integral 
in the exponent is performed over the whole classically forbidden domain, i.e. between the 
inner turning points of the orbits in the two wells. Preparing the initial state as one of the 
quasi-modes (i.e., as the even or odd superposition of the symmetric and the antisymmetric 
eigenstate) , the system will undergo Rabi oscillations between the wells with the frequency 
AE/h. The "tunneling rate" of this system is therefore given by the splitting (18. ip . Keeping 
all classical parameters flxed, it decreases exponentially with 1/h, and, in that sense, one 
can say that tunneling "vanishes" in the classical limit. 



8.1.2 Chaos-assisted tunneling 

The approach presented in the previous section can be generalized to multidimensional, even 
non-separable systems, as long as their classical dynamics is still integrable [S]. It breaks 
down, however, as soon as a non-integrable perturbation is added to the system, e.g. if the 
one- dimensional double-well potential is exposed to a driving that is periodic in time (with 
period r, say). In that case, the classical phase space of the system generally becomes a 
mixture of both regular and chaotic structures. 

As visualized by the stroboscopic Poincare section — which is obtained by retaining the 
phase space coordinates at every integer multiple of the driving period r — the phase space 
typically displays two prominent regions of regular motion, corresponding to the weakly per- 
turbed dynamics within the two wells, and a small (or, for stronger perturbations, large) 
layer of chaotic dynamics that separates the two regular islands from each other. Numerical 
calculations of model systems in the early nineties [HI HD] have shown that the tunnel split- 
tings in such mixed systems generally become strongly enhanced compared to the integrable 
limit. Moreover, they do no longer follow a smooth exponential scaling with 1/^ as expressed 
by Eq. (18. ip . but display huge, quasi-erratic fluctuations when h or any other parameter of 
the system vary [9l [10] . 

These phenomena are traced back to the speciflc role that chaotic states play in such 
systems [H], [121 [THl [M] . contrast to the integrable case, the tunnel doublets of the localized 
quasi-modes are, in a mixed regular-chaotic system, no longer isolated in the spectrum, but 
resonantly interact with states that are associated with the chaotic part of phase space. Due 
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to their delocalized nature in phase spac^, such chaotic states typically exhibit a significant 
overlap with the boundary regions of both regular wells. Therefore, they may provide an 
efficient coupling mechanism between the quasi-modes - which becomes particularly effective 
whenever one of the chaotic levels is shifted exactly on resonance with the tunnel doublet. 
As illustrated in Fig. 18. this coupling mechanism generally enhances the tunneling rate, 
but may also lead to a complete suppression thereof, arising at specific values of h or other 
parameters [T6] . 

We point out that this type of resonant tunneling does not necessarily require the presence 
of classical chaos and may appear also in integrable systems, for i nstance in a one-dimensional 
symmetric triple-well potential. This is illustrated in Fig. 18.21 which shows the scaling of 
level splittings associated with the two lateral wells as a function of 1 / h for such a triple- well 
potential. On top of an exponential decrease according to Eq. (18.11) . the splittings display 
strong spikes occuring whenever the energy of a state localized in the central well becomes 
quasi-degenerate with the energies of the states in the two lateral wells. 



e 




Figure 8.1: The two elementary scenarios of huge enhancement (case b) or cancellation (case 
a) of the tunneling splitting between symmetric (+) and anti-symmetric (— ) states can be 
easily understood from the resonant crossing of a third level (here a symmetric one) and 
the corresponding level repulsion between states of the same symmetry class. The external 
parameter that triggers the fluctuation of tunneling can be of quantum origin (effective h) 
or classical. 

For mixed regular-chaotic systems, the validity of this "chaos-assisted" tunneling picture 
was essentially confirmed by successfully modeling the chaotic part of the quantum dynamics 
with a random matrix from the Gaussian orthogonal ensemble (GOE) [HI |T2l |T7]. Using 
the fact that the coupling coefficients between the regular states and the chaotic domain are 
small, this random matrix ansatz yields a truncated Cauchy distribution for the probability 
density to obtain a level splitting of the size AE. Such a distribution is indeed encountered 
in the exact quantum splittings, which was demonstrated for the two-dimensional quartic 
oscillator [T7] as well as, later on, for the driven pendulum Hamiltonian that describes 



^Wc assume here that the effects of dynamical locahzation observed in Ref. 1151 remain irrelevant. 
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Figure 8.2: Resonant tunneling at work for a one dimensional triple well potential (here 
V{x) = {x^ — a^)^(x^ — fe^) with a = 1.75 and b = .5). When 1/h is increased, tunneling 
between the two lateral wells is enhanced by several orders of magnitude each time a level 
in the central well crosses down the doublet. 



the tunneling process of cold atoms in periodically modulated optical lattices [U [18]. A 
quantitative prediction of the average tunneling rate, however, was not possible in the above- 
mentioned theoretical works. As we shall argue later on, this average tunneling rate is directly 
connected to the coupling matrix element between the regular and the chaotic states, and 
the strength of this matrix element was unknown and introduced in an ad-hoc way. 

A first step towards this latter problem was undertaken by Podolskiy and Narimanov [T9] 
who proposed an explicit semiclassical expression for the mean tunneling rate in a mixed 
system by assuming a perfectly clean, harmonic-oscillator like dynamics within the regular 
island and a structureless chaotic sea outside the outermost invariant torus of the island. This 
expression turned out to be successful for the reproduction of the level splittings between 
near-degenerate optical modes that are associated with a pair of symmetric regular islands 
in a non-integrable micro-cavity [19] (see also Ref. |20J ) . The application to dynamical 
tunneling in periodically modulated optical lattices j^, for which splittings between the 
left- and the right-moving stable eigenmodes were calculated in Ref. [1], seems convincing 
for low and moderate values of l//i, but reveals deviations deeper in the semiclassical regime 
where plateau structures arise in the tunneling rates. Further, and more severe, deviations 
were encountered in the application of this approach to tunneling processes in other model 
systems |21j . 
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Backer, Ketzmerick, Lock, and coworkers [221 123] recently undertook the effort to de- 
rive more rigorously the regular-to-chaotic coupling rate governing chaos-assisted tunneling. 
Their approach is based on the construction of an integrable approximation for the noninte- 
grable system, designed to accurately describe the motion within the regular islands under 
consideration. The coupling rate to the chaotic domain is then determined through the com- 
putation of matrix elements of the Hamiltonian of the system within the eigenbasis of this 
integrable approximation |22] . This results in a smooth exponential-like decay of the aver- 
age tunneling rate with 1/h, which was indeed found to be in very good agreement with the 
exact tunnehng rates for quantum maps and billiards [221 123] . Those systems, however, were 
designed such as to yield a "clean" mixed regular-chaotic phase space, containing a regular 
island and a chaotic region which both do not exhibit appreciable substructures [22l [23]. 



8.1.3 The role of nonlinear resonances 

In more generic systems, such as the quantum kicked rotor or the driven pendulum [1], 
however, even the "average" tunneling rates do not exhibit a smooth monotonous behaviour 
with 1/h, but display peaks and plateau structures that cannot be accounted for by the 
above approaches. To understand the origin of such plateaus, it is instructive to step back 
to the conceptually simpler case of nearly integrable dynamics, where the perturbation from 
the integrable Hamiltonian is sufficiently small such that macroscopically large chaotic layers 
are not yet developed in the Poincare surface of section. In such systems, the main classical 
phase space features due to the perturbation consist in chain-like substructures that surround 
stable periodic orbits or equilibrium points of the classical motion. Those substructures 
come from nonlinear resonances between the internal degrees of freedom of the system or, 
for driven systems, between the external driving and the unperturbed oscillations around the 
central orbit. In a similar way as for the quantum pendulum Hamiltonian, such resonances 
induce additional tunneling paths in the phase space, which lead to couplings between states 
that are located near the same stable orbit [2H [25] . 

The relevance of this effect for the near-integrable tunneling process between two symmetry- 
related wells was first pointed out by Bonci et al. [26] who argued that such resonances may 
lead to a strong enhancement of the tunneling rate, due to couplings between lowly and 
highly excited states within the well which are permitted by near-degeneracies in the spec- 
trum In Refs. [271 [28] , a quantitative semiclassical theory of near-integrable tunneling was 
formulated on the basis of this principal mechanism. This theory allows one to reproduce the 
exact quantum splittings from purely classical quantities and takes into account high-order 
effects such as the coupling via a sequence of different resonance chains [271 EH] • More re- 
cent studies by Keshavamurthy on classically forbidden coupling processes in model systems 
that mimic the dynamics of simple molecules confirm that the resonance-assisted tunneling 
scenario prevails not only in one-dimensional systems that are subject to a periodic driving 
(such as the kicked Harper model studied in Re f. [271 [28]). but also in autonomous systems 
with two and even three degrees of freedom [291 [30] . 

In Refs. [31] [32] [33] [M] resonance-assisted couplings were incorporated in an approximate 
manner into the framework of chaos-assisted tunneling in order to provide a quantitative the- 
ory for the regular-to-chaotic coupling rate. In this context, it is assumed that the dominant 
coupling between regular states within and chaotic states outside the island is provided by 
the presence of a nonlinear resonance within the island. A straightforward implementation 
of this idea yields good agreement with the exact tunneling rates as far as their average 
decay with 1/^ in the deep semiclassical limit is concerned. Moreover, individual plateaus 
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and peak structures could be traced back to the influence of specific nonlinear resonances, 
not only for double- well- type tunneling in closed or periodic systems [311 |32l |33], but also 
for tunneling- induced decay in open systems [31]. However, the predictive power of this 
method was still rather limited, insofar as individual tunneling rates at given system pa- 
rameters could be over- or underestimated by many orders of magnitude. In particular, 
resonance-assisted tunneling seemed inapplicable in the "quantum" limit of large h, where 
direct regular-to chaotic tunneling proved successful 



A major advance in this context was achieved by improving the semiclassical evaluation 
of resonance-induced coupling processes in mixed systems, and by combining it with "direct" 
regular-to-chaotic tunneling [35]. This combination resulted, for the first time, in a semi- 
classical prediction of tunneling rates in generic mixed regular-chaotic systems that can be 
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compared with the exact quantum rates on the level of individual peak structures [35 
confirms the expectation that nonlinear resonances do indeed form the "backbone 
non-monotonous substructures in tunneling rates. It furthermore suggests that those rates 
could, also in systems with more degrees of freedom, possibly be estimated in a quantita- 
tively satisfactory manner via simple classical computations, based on the most prominent 
nonlinear resonances that are manifested within the regular island. 

It is in the spirit of this latter expectation that this contribution has been written. Our 
aim is not to formulate a formal semiclassical theory of tunneling in mixed systems, which 
still represents an open problem that would rather have to be solved on the basis of complex 
classical orbits [361 EZl [38] • Instead, we want to provide a simple, easy-to-implement, yet 
effective prescription how to compute the rates and time scales associated with tunneling 
processes solely on the basis of the classical dynamics of the system, without performing 
any diagonalization (not even any application) of the quantum Hamiltonian or of the time 
evolution operator. This prescription is based on chaos- and resonance-assisted tunneling in 
its improved form [35] . The main part of this contribution is therefore devoted to a detailed 
description of r eson ance -ass isted tunneling and its combination w ith c haos-assisted tunneling 
in the sections 18.21 and 18. 3[ respectively. We present in section 18.41 the application of this 
method to tunneling processes in the quantum kicked rotor, and discuss possible limitations 
and future prospects in the conclusion in section 18.51 



8.2 Theory of resonance-assisted tunneling 



8.2.1 Secular perturbation theory 



For our study, we restrict ourselves to systems with one degree of freedom that evolve under 
a periodically time-dependent Hamiltonian H{p,q,t) = H{p,q,t + t). We suppose that, 
for a suitable choice of parameters, the classical phase space of H is mixed regular-chaotic 
and exhibits two symmetry-related regular islands that are embedded within the chaotic sea. 
This phase space structure is most conveniently visualized by a stroboscopic Poincare section, 
where p and q are plotted at the times t = nT{n G Z). Such a Poincare section typically 
reveals the presence of chain-like substructures within the regular islands, which arise due 
to nonlinear resonances between the external driving and the internal oscillation around the 
island's center. Before considering the general situation for which many resonances may 
come into play in the tunneling process, we start with the simpler case where the two islands 
exhibit a prominent r:s resonance, i.e., a nonlinear resonance where s internal oscillation 
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periods match r driving periods and r sub-islands are visible in the stroboscopic section. 

The classical motion in the vicinity of the r:s resonance is approximately integrated by 
secular perturbation theory [32] (see also Ref. For this purpose, we formally introduce 

a time- independent Hamiltonian Hq{p, q) that approximately reproduces the regular motion 
in the islands and preserves the discrete symmetry of H. In some circumstances, as for 
instance if H is in the nearly integrable regime, Ho{p,q) can be explicitly computed within 
some approximation scheme (using for instance the Lie transformation method [39]). We 
stress though that this will not always be necessary. Assuming the existence of such a Hq, the 
phase space generated by this integrable Hamiltonian consequently exhibits two symmetric 
wells that are separated by a dynamical barrier and "embed" the two islands of H. In terms 
of the action-angle variables (/, 6) describing the dynamics within each of the wells, the total 
Hamiltonian can be written as 

H{I,e,t) = Ho{I) + V{I,e,t) (8.2) 
where V would represent a weak perturbation in the center of the island 

The nonlinear r:s resonance occurs at the action variable Ir.s that satisfies the condition 

rflr:s = su (8.3) 

with u = 2'it/t and 

a.. = ^ . (8.4) 



dl 



I = Ir: 



We now perform a canonical transformation to the frame that corotates with this resonance. 
This is done by leaving I invariant and modifying 6 according to 

0^^ = 0-Q^.j_ (8.5) 

This time-dependent shift is accompanied by the transformation H = H — fi^.^J in 
order to ensure that the new corotating angle variable i} is conjugate to /. The motion of I 
and 'd is therefore described by the new Hamiltonian 

■H{I,^,t)=no{I) + V{I,^,t) (8.6) 

with 

Hoil) = Ho{I) - nr:sl , (8.7) 

V{I,^,t) = V{I,^ + nr:st,t) . (8.8) 



The expansion of Ho in powers of / — Ir:s yields 

noil) ^ + ^^^f +0[{I- Ir:s?] (8.9) 

with a constant H^^^ = Ho{Ir:s) — ^Ir-.slr-.s and a quadratic term that is characterized by 
the effective "mass" parameter rrir-.s = [(PHo/ dP{Ir:s)\~^ ■ Hence, dT-Lo/dl is comparatively 
small for I ~ 7^.^, which implies that the co-rotating angle varies slowly in time near 

^In order not to overload the notation, we use the same symbol H for the Hamiltonian in the original phase-space variables 
(p, q) and in the aetion-angle variables (/, 8). 
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the resonance. This justifies the apphcation of adiabatic perturbation theory [39|, which 
effectively amounts, in first order, to replacing V(/, t) by its time average over r periods 
of the driving (using the fact that V is periodic in t with the period rr) 0. We therefore 
obtain, after this transformation, the time-independent Hamiltonian 

H{I,^)=Ho{I) + V{I,^) (8.10) 

with 

V(/,^9) = — / V{I,'d,t)dt. (8.11) 
rr Jo 

By expanding V{I, 6, t) in a Fourier series in both 6 and t, i.e. 

oo 

V{I,e,t)= ^ V^,„,(/)e*'V™"* (8.12) 

i,m=— oo 

with Vi^m{I) = [V-i-m{I)]*, one can straightforwardly derive 

oo 

V(/, ^) = Vo,o{I) + Yl 2^fc(^) ^os{kr^ + 0,) (8.13) 

k=l 

defining 

Vk{I)e'^''=Vrk,-sk{I), (8.14) 
i.e., the resulting time-independent perturbation term is (27r/r)-periodic in i). 

For the sake of clarity, we start discussing the resulting effective Hamiltonian neglecting 
the action dependence of the Fourier coefficients of V{I,'&). We stress that this dependence 
can be implemented in a rel atively straightforward way using Birkhoff-Gustavson normal- 
form coordinates (cf section [8.2.31 below): it is actually important to obtain a good quan- 
titative accuracy. For now, however, we replace Vk{I) by = Vfc(/ = Ir.s) in Eq. fl8.13p . 
Neglecting furthermore the term Vo,o(-^)) we obtain the effective integrable Hamiltonian 

oo 
k=l 

for the description of the classical dynamics in the vicinity of the resonance. We shall see in 
section [8.2.21 that the parameters of H^^s relevant to the tunneli ng pro cess can be extracted 
directly from the classical dynamics of H{t), which is making Eq. (18.15!) particularly valuable. 



8.2.2 The pendulum approximation 

The quantum implications due to the presence of this nonlinear resonance can be straight- 
forwardly inferred from the direct semiclassical quantization of Hj-^s, given by 

oo 

4es = ^o(/) - ^r J + cos(A;r^^ + 0^) . (8.16) 

k=l 

■'This step involves, strictly speaking, another time-dependent canonical transformation {I, i9) ^ (/, i9) which slightly modifies 
/ and i9 (see also Ref. [28]). 
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Here we introduce the action operator I = —itid/dd and assume anti-periodic boundary 
conditions in -d in order to properly account for the Maslov index in the original phase space 
|24] . In accordance with our assumption that the effect of the resonance is rather weak, 

we can now apply quantum perturbation theory to the Hamiltonian fl8.16p . treating the I- 

dependent "kinetic" terms as unperturbed part and the -i^-dependent series as perturbation. 
The unperturbed eigenstates are then given by the (anti-periodic) eigenfunctions {d\n) = 

(27r)~^/^ exp[«(n + l/2)'i?] {n > 0) of the action operator / with the eigenvalues 

4 = % + 1/2). (8.17) 



As is straightforwardly evaluated, the presence of the perturbation induces couplings 
between the states \n) and \n + kr) with the matrix elements 

(n + kr\H,^,\n) = Vke''^" (8.18) 

for strictly positive integer k. As a consequence, the "true" eigenstates \ipn) of Hres contain 
admixtures from unperturbed modes \n') that satisfy the selection rule \n' — n\ = kr with 
integer k. They are approximated by the expression 

I / \ I \ , {n + kr\H,es\n) , , \ , 
\w„) = \n) + y \n + kr) + 

' ^ ' ^ E^-Er^+kr + kshw^ ' 

^ {n + kr\H,,,\n + k'r) {n + k'r\H,,,\n) , , _ _ _ /g^^g) 

En- En+kr + kstUxl E„ - E^+k'r + k' stllO 

where En = Ho{In) denote the unperturbed eigenenergies of Hq and the resonance condition 
(18. 3p is used. The summations in Eq. (18.191) are generally finite due to the finiteness of the 
phase space area covered by the regular region. 

W ithin the quadratic approximation of Hq{I) around Ir s, we obtain from Eqs. (18. 7p and 

(ESI) 

En - Ho{Ir,s) + ^r:s{In " Ir:s) + 77^(4 " ^r-sf • (8.20) 



2mr- 



r:s 



This results in the energy differences 

En - En+kr + kshw ~ 2m~^^'^ " In+kr){In + h.+kr " '^Ir.s) ■ (8.21) 

From this expression, we see that the admixture between \n) and \n') becomes particularly 
strong if the r:s resonance is symmetrically located between the two tori that are associated 
with the actions and — i.e., if -|- /„/ ^ ^I^.s- The presence of a significant nonlinear 
resonance within a region of regular motion provides therefore an efficient mechanism to 
couple the local "ground state" — i.e, the state that is semiclassically localized in the center 
of that region (with action variable /q < Ir.s) — to a highly excited state (with action 
variable Ikr > Ins)- 

It is instructive to realize that the Fourier coefficients Vk of the perturbation operator 
decrease rather rapidly with increasing k. Indeed, one can derive under quite general cir- 
cumstances the asymptotic scaling law 

Vk ~ {kr)Woexp[-krQr;stim{Ir:s)] (8.22) 



10 



CHAPTER 8. RESONANCE-ASSISTED TUNNELING 



for large k, which is based on the presence of singularities of the complexified tori of the 
integrable approximation Hq{I) (see Eq. (66) in Ref. [2H])- Here tim(-^) denotes the imaginary 
time that elapses from the (real) torus with action / to the nearest singularity in complex 
phase space, 7 corresponds to the degree of the singularity, and Vq contains information about 
the correspondi ng res idue near the singularity as well as the strength of the perturbation. 
The expression fl8.22p is of little practical relevance as far as the concrete determination of 
the coefficients Vk is concerned. It permits, however, to estimate the relative importance 
of different perturbative pathways connecting the states \n) and \n + kr) in Eq. fl8.19p . 
Comparing e.g. the amplitude A2 associated with a single step from \n) to \n + 2r) via V2 
and t he am plitu de Ai associated with two steps from \n) to \n + 2r) via Vi, we obtain from 
Eqs. flOB and fl8:22ll the ratio 

mr-.sVo 

under the assumption that the resonance is symmetrically located in between the correspond- 
ing two tori (in which case we would have In+r — h-.s)- Since Vq can be assumed to be finite 
in mixed regular-chaotic systems, we infer that the second-order process via the stronger 
coefficient Vi will more dominantly contribute to the coupling between |n) and \n + 2r) in 
the semiclassical limit ^ — 0. 

A similar result is obtained from a comparison of the one-step process via Vk with the 
fc-step process via Vi, where we again find that the latter more dominantly contributes to 
the coupling between \n) and \n + kr) in the limit — ?■ 0. We therefore conclude that in 
mixed regular-chaotic systems the semiclassical tunneling process is adequately described by 
the lowest non- vanishing term of the sum over the Vk contributions, which in general is given 
by Vi cos(r'i9-|-0i)0. Neglecting all higher Fourier components Vk with k > 1 and making the 
quadratic approximation of Hq around / = Ir-s-, we finally obtain an effective pendulum-like 
Hamiltonian ^ 

i7res(/, ^) ^ ^^'^'"^ + 2Vr:s COs{rd + 0i) (8.24) 

with Vr:s = Vi^. 

This simple form of the effective Hamiltonian allows us to determine the parameters Ir.s, 
rrir-.s and Vr-.s from the Poincare map of the classical dynamics, without explicitly using the 
transformation to the action- angle variables of Hq. To this end, we numerically calculate 
the monodromy matrix M^.^ = d{pf, qf)/d{pi, qi) of a stable periodic point of the resonance 
(which involves r iterations of the stroboscopic map) as well as the phase space areas S^^ 
and S~.g that are enclosed by the outer and inner separatrices of the resonance, respectively 
(see also Fig. 18.31) . Using the fact that the trace of Mr-.s as well as the phase space areas S^^ 
remain invariant under the canonical transformation to {!,'&), we infer 



Ir:s = — (5+ +5-J, (8.25) 

V^mr-.sVr-.s = ^iS^.,-S-,), (8.26) 
16 

/2K^ = ^arccos(trM^.,/2) (8.27) 
V rrir-.s r'^T 



^Exceptions from this general rule typically arise in the presence of discrete symmetries that, e.g., forbid the formation of 
resonance chains with an odd number of sub-islands and therefore lead to Vi = for an r:s resonance with an odd r. 
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from the integration of the dynamics generated by i/res [10] • Eqs. fl8.25p - fl8.27p make it 
possible to derive the final expressions for the tunneling rates directly from the dynamics of 
H{t), without explicitly having to construct the integrable approximation Hq and making 
the Fourier analysis of V{p, q, t) = H{p, q, t) — HQ{p, q). As this construction of the integrable 
approximation may turn out to be highly non-trivial in the mixed regime, avoiding this step 
is actually an essential ingredient to make the approac h we are following practical. We note 
though that improving the quadratic approximation fl8.24p for Hq is sometimes necessary, 
but this does not present any fundamental difficulty. 




Figure 8.3: Classical phase space of the kicked rotor Hamiltonian aX K = 3.5 showing a 
regular island with an embedded 6:2 resonance. The phase space is plotted in the original 
(p, q) coordinates (u pper l eft panel), in approximate normal-form coordinates (P, Q) defined 
by Eqs. fl8.28p and f l8.29p (upper right panel), and in approximate action-angle variables 
{I^'d) (lower panel). The thick solid and dashed lines represent the "outer" and "inner" 
separatrix of the resonance, respectively. 
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8.2.3 Action dependence of the coupling coefficients 



Up to now, and in our previous publications [271 ESI EB [321 133] , we completely neglected the 
action dependence of the coupling coefficients Vk{I). This approximation should be justified 
in the semiclassical limit of extremely small h, where resonance-assisted tunneling generally 
involves multiple coupling processes [28] and transitions across individual resonance chains 
are therefore expected to take place in their immediate vicinity in action space. For finite 
h, however, the replacement Vfc(J) i— ?■ Vk{Ir:s), permitting the direct quantization in action- 
angle space, is, in general, not sufficient to obtain an accurate reproduction of the quantum 
tunneling rates. We show now how this can be improved. 

To this end, we make the general assumption that the classical Hamiltonian H{p,q,t) of 
our system is analytic in p and q in the vicinity of the regular islands under consideration. 
It is then possible to define an analytical canonical transformation from {p, q) to Birkhoff- 
Gustavson normal-form coordinates (P, Q) [HI |12] that satisfy 



P = -V2Jsin^, (8.28) 

Q = V2Icos9 (8.29) 

and that can be represented in power series in p and q. The "unperturbed" integrable 
Hamiltonian Hq therefore depends only on / = (P^ + Q'^)/2. 



Writing 

V V2I 

for positive /, we obtain, from Eq. fl8.12p . the series 



^.30) 



oo oo 

vii,e,t)= + + + [ 

m=— oo I. 1=1 V 2/ J 



(8.31) 

for the perturbation. Using the fact that V{I, 9, t) is analytic in P and Q, we infer that 
Vi^rn{I) must scale at least as /'/^. By virtue of fl8.14p . this implies the scaling Vk{I) oc P^l"^ 
for the Fourier coefficients of the time-independent perturbation term that is associated with 
the r:s resonance. Making the ansatz Vfc(/) = P^^'^Vk (and neglecting the residual action 
dependence of Vk), we rewrite Eq. f l8.13p a^ 

Vk 



V(/, ^) = Vo,o{I) + E ^ [(^ - ^^)''^'^' + + ^Pf^'e-'^'] . (8.32) 



k=l 



Each term in the sum is given by the well-known Birkhoff normal form generically describing 
a r > 3 resonance when the bifurcation of the stable periodic orbit is controlled by one single 
parameter (see, e.g., Eq. (4.70) in Chapter 4 of Ref. [13] or Eqs. (3.3.17) and (3.3.18) in 
Ref. [S] for another simple derivation of the action dependence of Vfc(/)). The term Vo,o(-^) 
is neglected in the following as it does not lead to any coupling between different unperturbed 
eigenstates in the quantum system. 



^This involves, strictly speaking, another canonical transformation of P and Q to the frame that co-rotates with the resonance. 
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One comment is in order here. The small parameter in th e adiabatic approximation (18.101) 
is the difference (I—Ir-.s), while in the derivation of Eq. fl8.32p we have neglected higher powers 
of / and thus assumed the action I itself to be small (strictly speaking we should work 
near one bifurcation only). We thus mix a development near the resonant torus with one 
near the center of the island. This may eventually become problematic if (i) the resonance 
chain is located far away from the center of the island, in which case the associated coupling 
strength may contain a nonnegligible relative error when being computed via the assumption 
Vfc(/) = F^/'^Vk with constant -Ufc, and if (ii) that coupling strength happens to appear rather 
often in the main perturbative chain that connects the quasimodes of the island to the 
chaotic sea, which generally would be the case for low-order resonances with relatively small 
r3 Otherwise, we expect that this inconsistency in the definition of the regimes of validity 
of our perturbative approach does not lead to a significant impact on the numerical values of 
the semiclassical tunneling rates, which indeed seems to be confirmed by numerical evidence 
to be discussed below. 

This being said, the quantization of the resulting classical Hamiltonian can now be carried 
out in terms of the "harmonic oscillator" variables P and Q and amounts to introducing the 
standard ladder operators a and according to 

^ {Q + iP), (8.33) 



2h 

at = ^(g-^P). (8.34) 



This yields the quantum Hamiltonian 

oo 

ifres = H^{i) - ^r J + Vkti""-'^ [a'^^c"^'^^ + (at)'=^e^'^^] (8.35) 



fc=i 



with I = /i(ata + l/2). As for Eq. (18.161) . perturbative couplings are introduced only between 
unperturbed eigenstates \n) and \n') that exhibit the selection rule \n' — n\ = kr wi th integer 
k. The associated coupling matrix elements are, however, different from Eq. (I8.18P and read 



{n + kr\H,M = C,y^'V<^'=ji!l±M: 

V n\ 



V,(I..,.)e*[^\"\[^^ (8,36) 



for strictly positive k. Close the resonance, i.e. more formally taking the semiclassical 
limit 77, — )■ oo keeping k and 5 = (Ir-.s/f^ — n) fixed, and making use of the Stirling formula 

n\ ~ A/2'/rn(n/e)", Eq. (I8.36P reduces to Vk^Ir-.s)^^'^'' ■ The difference becomes, on the other 
hand, particularly pronounced if the r:s resonance is, in phase space, rather asymmetrically 
located in between the invariant tori that correspond to t he sta tes \n) and \n + kr) — i.e., 
if Ir:s is rather close to In or to In+kr- In that case, Eq. (I8.18P may, respectively, strongly 
over- or underestimate the coupling strength between these states. 



^ Corrections to the form 118.321 1 should also arise in the presence of prominent secondary resonances, which occur when 
primary resonances start to overlap and create chains of sub-islands nested inside the primary island chains. 
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8.2.4 Multi-resonance processes 



Up to this point, we considered the couphngs between quasi- modes generated by a given 
resonance. In general, however, several of them may play a role for the coupling to the 
chaotic sea, giving rise to multi-resonance transitions across subsequent resonance chains 
in phase space [271 EH]- As was argued in the context of near-integrable systems such 
multi-resonance processes are indeed expected to dominate over couplings involving only one 
single resonance in the deep semiclassical limit ^ — 0. 

The description of the coupling process across several consecutive resonances requires 
a generalization of Eq. (18.1 9p describing the modified eigenstate due to resonance-induced 
admixtures. We restrict ourselves, for this purpose, to including only the first-order matrix 

elem ents {n + r|^res^''|?^) for each resonance [i.e., only the matrix elements with k = 1 in 
Eq. (I8.36P ]. For the sake of clarity, we furthermore consider the particular case of coupling 
processes that start in the lowest locally quantized eigenmode with node number n = 
(the generalization to initial n ^ being straightforward). The prescription we use is to 
consider that, although the approximation (IS.lSp is valid for only one resonance at a time, 
it is possible to sum the contributions obtained from different resonances. Considering a 
sequence of consecutive r:s, r':s', r":s" . . .resonances, we obtain in this way 



i«>+E n 

fc>0 \/=i 



(/r|i74T)|(/-l)r) 



En — E, 



Ir 



ishjjj 



X 



k' 



iH+E n 



{kr + l'r'\H'^i"'\kr + {V - 1)/) 
Eq — Ekr+i'r' + {ks + V s')fyjj 

k" 



k'>0 \l 



X 



X 



\kr 



, , ^ I TT {kr + k'r' + l"r"\H^C-'"\kr + k'r' + (/" - l)r") 
k'r') + E n 



k">0 \l"=l 



Eq — Ekr+k'r'+l"r" + {ks + I's' + l"s")huj 



X 



X 



\kr + k'r' + k"r") + ^ 



fc"'>0 



(8.37) 



for the modified "ground state" within the island. Given an excit ed qua si-mode n far from 
the interior of the island, the overlap {nlipo) obtained from Eq. (18.371) will in most cases 
be exponentially dominated by one or a few contributions. There is no systematic way to 
identify them a priori, although some guiding principle can be used in this respect 



Quite naturally, for instance, low-order resonances, with comparatively small r and s, will, 
in general, give larger contribution than high-order resonances with comparable winding 
numbers s/r but larger r and s, due to the strong differences in the sizes of the mean 



coupling matrix elements Vr-.s [see, e.g., Eq. 



In the same way, sequences of couplings 



involving small denominators, i.e. energy differences like Eq. (I8.2ip that are close to zero, and 
thus intermediate steps symmetrically located on each side of a resonance, will tend to give 
larger contributions. In the small h limit this will tend to favor multi-resonance processes. 
Conversely, for intermediate values of h (in terms of the area of the regular region) the main 
contributions can be obtained from the lowest-order resonances. With few exceptions — 
especially concerning low-order resonances that are located close to the center of the island, 
thereby leading to relatively large energy denominators and small admixtures — this rule 
is generally observed for the semiclassical calculation of the eigenphase splittings we shall 
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8.3 The combination with chaos-assisted tunnehng 



We now discuss the implication of such nonhnear resonances on the tunnehng process between 
the two symmetry-related regular islands under consideration. In the quantum system, these 
islands support (for not too large values of h) locally quantized eigenstates or "quasimodes" 
with different node numbers n, which, due to the symmetry, have the same eigenvalues 
in both islands. In our case of a periodically driven system, these eigenvalues can be the 

eigenphases ipn of the unitary time evolution (Floquet) operator U over one period r of the 
driving, or, alternatively, the quasienergies En such that = —Enr/h (modulo 27r). 

The presence of a small (tunneling-induced) coupling between the islands lifts the de- 
generacy of the eigenvalues and yields the symmetric and antisymmetric linear combination 
of the quasimodes in the two islands as "true" eigenstates of the system. A nonvanishing 
splitting A<y9„ = \lp^ — lp~\ consequently arises between the eigenphases ip^ of the symmetric 
and the antisymmetric state, which is related to the splitting AE'^ = \E^ — E~\ of the 
quasi-energies E^ through A<y9„ = rAEn/h. 



8.3.1 Resonance-assisted tunneling in near-integrable systems 

We start by considering a system in the nearly integrable regime. In that case, we can 
assume the presence of a (global) integrable Hamiltonian Ho{p, q) that describes the dynamics 
in the entire phase space to a very good approximatioiiH- The energy splittings for the 
corresponding quantum Hamiltonian Hq = HQ{p,q) can be semiclassically calculated via an 
analytic continuation of the invariant tori to the complex domain [8] . This generally yields 
the sphttings 



(up to a numerical factor of order one) where Qn is the classical oscillation frequency associ- 
ated with the nth quantized torus and 0"„ denotes the imaginary part of the action integral 
along the complex path that joins the two symmetry-related tori. 

The main effect of nonlinear resonances in the non-integrable system is, as was discussed 
in the previous subsections, to induce perturbative couplings between quasimodes of different 
excitation within the regular islands. For the nearly integrable systems this can already lead 
to a substantial enhancement of the splittings AEn as compared to Eq. fIS.SSp [27tl28]. As can 
be derived within quantum perturbation theory, the presence of a prominent r:s resonance 
modifies the splitting of the local "ground state" in the island (i.e., the state with vanishing 



^Formally, this Hamiltonian is not identical with the unperturbed approximation Ho{I) introduced in section [8.2.11 as the 
definition of the latter is restricted to one well only. It is obvious, however, that Ho{I) can be determined from Ho(p,q), e.g. 
by means of the Lie transformation method 1391 . 
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node number n = 0) according to 

Av.o = Av.r + i:i47^PAv.i? (8.39) 

k=l 

(using Aifn^ ^ ^'^o'^ for n > 0), where A^^^'^^ = (krltpo) denotes the admixture of the 
(fcr) th exc ited unperturbed co mpon ent \kr) to t he pe rturbed ground state lipo) according to 
Eq. fl8.19p [possibly using Eq. fl8.36p instead of f l8.18p ]. The maximal number kc of coupled 
states is provided by the finite size of the island according to 



k 



c 



1 / area of the island 1 
r V ^ 2 



(8.40) 



(r-.s) 
kr 



where the bracket stands for the integer part. The rapid decrease of the amplitudes A 

with k is compensated by an exponential increase of the unperturbed splittings A(^^°^ , arising 
from the fact that the tunnel action (j„ in Eq. f l8.38p generally decreases with increasing n. 
The maximal contribution to the modified ground state splitting is generally provided by the 
state \kr) for which I^r + Io — '^Ir-.s — i-e., which in phase space is most closely located to the 
torus that lies symmetrically on the opposite side of the resona nce c hain. This contribution is 
particularly enhanced by a small energy denominator [see Eq. fl8.2ip ] and typically dominates 
the sum in Eq. fl8.39p . 

As one goes further in the semi classical ^ — >■ limit, a multi-resonance process is usually 
the dominant one. Neglecting interference terms between different coupling pathways that 
connect the ground state with a given excited mode \n) (which is justified due to the fact 
that the amplitudes associated with those coupling pathways are, in general, much different 
from each other in size), we obtain from Eq. fl8.37p an expression of the form 

A^o = A^r + Y: + E E I VKU'PA^£U' + • • • (8-41) 

k k k' 

with the coupling amplitudes 



(/r|g,^eT1(/-l)r) 

Eq — Eir + Ishu 



- 11 R. ^/,.,„ (8-42) 



1=1 



Ar':s') _ ij {kr + Vr'\m:.-"\kr + {l'-iy) 

" Eo~E,r+i'r' + {ks + l's')hu; ^'-^'^ 



•^kr+k'r',kr+k'r'+k"r" 

for the eigenphase splitting. 



8.3.2 Coupling with the chaotic sea 



Turning now to the mixed regular-chaotic case, the integrable Hamiltonian Hq{I) provides 
a good approximation only near the center of the regular island under consideration, and 
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invariant tori exist only up to a maximum action variable Ic corresponding to the outermost 
boundary of the regular island in phase space. Beyond this outermost invariant torus, mul- 
tiple overlapping resonances provide various couplings and pathways such that unperturbed 
states in this regime can be assumed to be strongly connected to each other. Under such 
circumstances, it is natural to divide the Hilbert space into two parts, integrable and chaotic, 
associated respectively with the phase space regions within and outside the regular island. 

For each symmetry class ± of the problem, let us introduce an effective Hamiltonian H^ 
modeling the tunneling process. Let us furthermore denote Preg and Pch the (orthogonal) pro- 
jectors onto the regular and chaotic Hilbert spaces. The diagonal blocks H^^ = Pj-cg-f/^-Prog 
and H^ = Pch -ff^ Pch receive a natural interpretation: within H^^, on the one hand, the 

dynamics is exactly the same as in the nearly integrable regime above; H^, on the other 
hand, is best modeled in a statistical manner by the introduction of random matrix ensem- 
bles. The only remaining delicate point is thus to connect the two, namely to model the 

off-diagonal block P,.egH^Pch- We stress that there is not yet a real consensus on the best 
way how to do this, although various approaches give good quantitative accuracy. 

To state more clearly the problem, let us consider a regular state \n) with quasi-energy 
E^ close to the regular-chaos boundary. (Note that "close" here means that no resonance 
within the island can connect \n) to a state \n') within the island with n' > n. This notion of 
"closeness" to the boundary is therefore /i-dependent.) The resonance assisted mechanism 
will connect any quasi-mode deep inside the island to such a state at the edge of the island. 
But to complete the description of the chaos-assisted tunneling process it is necessary to 
compute the variance of the random matrix elements Vni between \n) and the eigenstates 

of (the variance is independent of i if H"^ is modeled by the Gaussian orthogonal 
or unitary ensemble). 

One possible approach to compute this quantity is the fictitious integrable system ap- 
proach that was proposed by Backer et al. [22] • This method relies on the fact that, for the 

effective Hamiltonian i^eff, the "direct" transition rate from a regular state |?t,) to the chaotic 
region is given using Fermi's golden rule by 



(see e.g. section 5.2.2 of Ref. [TT] for a discussion in the context of random matrix theory) 
where Ach denotes the mean spacing between eigenenergies within the chaotic block. As a 
consequence, one obtains in first order in ri;^//lAch, 



(see also the contribution of Backer, Ketzmerick, and Lock in this book). If one can explicitly 
construct a good integrable approximation ifreg = Ho{p,q) of the time- dependent dynamics 
(see Sec. I8.2.ip . this allows one, by quantum or semiclassical diagonalization, to determine 
the unperturbed eigenstates |?7,) within the regular island, and to construct the projectors Preg 

and Pch The "direct" regular-to-chaotic tunneling matrix elements of the nth quantized state 
within the island is then evaluated by a simple application of the quantum time evolution 

operator U over one period of the driving. 




h Ach 
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This approach can be quahfied as "semi-numerical" , as it requires to numerically perform 
the quantum evolution for one period of the map (although this can be done by analyti- 
cal and semiclassical techniques in some cases [23], see also the contribution of Backer et 

al. in this book). It strongly relies on the quality of the integrable approximation if^cg of 
the Hamiltonian. If the latter was really diagonal (which, as a matter of principle, cannot 
be achieved by means of classical perturbation theory, due to the appearance of nonlinear 
resonances), Eq. fl8.45p would represent an exact result (apart from the first-order approxi- 
mation in Tv'^/hAch which should not be a limitation in the tunneling regime). And indeed 
very good agreement between this prediction and numerically computed tunneling rates was 
found for quantum maps that were designed such as to yield a "clean" mixed regular-chaotic 
phase space, containing a regular island and a chaotic region which both do not exhibit 
appreciable substructures [22] , as well as for the mushroom billiard [23] . In more generic sit- 
uations, where nonlinear resonances are manifested within the regular island, this approach 
yields reliable predictions for the direct tunneling of regular states at the regular-chaos bor- 
der, and its combination with the resonance assisted mechanism described in the previous 
sections leads to good quantitative predictions for the tunneling rates for the states deep in 
the regular island [35] . We shall illustrate this on the example of the kicked rotor system in 
section 18.41 



8.3.3 Integrable semiclassical models for the regular-to-chaotic coupling 



For now, however, we shall discuss other possible approaches of purely semiclassical nature 
(i.e. not involving any numerical evolution nor diagonalization of the quantum system) 
to the calculation of the coupling parameter for a regular state 1^) at the edge of the 
regular-chaos boundary. 

One way to obtain an order of magnitude of is to consider the decay of the quasimode 
inside the regular island. For this purpose, let us assume that an integrable approximation 
Ho{I), valid up to a maximum action Jc corresponding to the chaos boundary, has be en ob- 
tained. Within the Birkhoff-Gustavson normal-form coordinates (P, Q) given by Eqs. (18.281) 
and (I8.29p . Hq appears as a function of the harmonic-oscillator Hamiltonian (P^ -|- Q^)/2 
and therefore has the same eigenstates 

{Q\n) = ^2)1/^ exp {-Q'/2n) UniQ/h'/') (8.46) 

^ ^ exp(- riPnmdo) for Q > Q„ , (8.47) 

where T-Ln are the Hermite polynomials and PniQ) = \/2In — Q^- The last equation cor- 
responds to the semiclassical asymptotics in the forbidden region on the right-hand side 
of the turning point Qi = \f2J^ with In = h{n -|- 1/2). In Q representation, the regular 
island extends up to Qc = \p2Tc-, at which point de state \n) has decayed to 4'n{Qc) — 
^^^^^ exp (-5(4, /J), where 

S{h,Ic)= \PnmdQ = VWc - In) - /nln ((^/c " 4 + Vh) / Vh) • (8.48) 
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This expression is suspiciously simple (as it depends only on 1^ and Jc, and on no other 
property of the system), and should not be taken too seriously as it is. 

Indeed, it should be borne in mind that f ? is related not so much to the value of the 
wavefunction at the regular-chaos boundary than to the transition rate r^_j^j,jj^Qg through 
Eq. (18.441) . which we can equal to the current flux Jn through this boundary for the regular 
state Using Hq{I) to compute this current leads to a zero result and it is therefore 
mandatory to use a better approximation of the non-integrable Hamiltonian H to obtain a 
meaningful answer. What complicates the evaluation of the transition rate from the edge 
of the regular to the chaotic domain is therefore that one needs to find an approximation 
describing both the regular and chaotic dynamics — unless one actually uses there the exact 
quantum dynamics as was done in Ref. p2]. Since the regular-chaos border is typically 
the place where approximation schemes tend to be difficult to control, this will rely on some 
assumption to be made for the chaotic regions, two possible choices of which we shall describe 
now. 

One scenario that has been considered amounts in some way to model the regular to chaos 
transition in the way depicted in Fig. 18.4b .: a kinetic-plus-potential Hamiltonian p^/2m+\^(g) 
where the island itself correspond to a potential well and the edge of the regular region to 
the place where the potential decreases abruptly. The picture one has in mind in that case 
is that escaping from the edge of the regular island to the chaotic sea is akin to the standard 
textbook barrier tunneling j35j. Using hanger's connection formula [16] within that model, 
the semiclassical wavefunction for the quasi-mode n inside the potential well can be extended 
under the potential barrier and, beyond this, into the region where motion at energy E'^ is 
again classically authorized. In the classically allowed region outside the well [q > q[.) the 
semiclassical wavefunction can be written as 



^^^^^ " V ^^'"P /%.(^)^^ + -/4) exp (-|) , (8.49) 
with the angular frequency of the torus E^, Pn{q) = A/2m[_E° — V{q)], and 

St= r\Pn{q)\dq (8.50) 

J qr 

the tunneling action^. The current of probability leaving the well is then given by 

J, ^ A Im[^* V^,] = ^ exp (-2|) , (8.51) 

from which is obtained through Eq. fl8.44p (identifying Jn with T'^^^^slos)- 

Although Eq. (18.511) is derived here for the particular case of a kinetic-plus-potential 
Hamiltonian, it applies more generally (up maybe to factors o f or der one) to any system 
with a phase space portrait that is similar to the one of Fig. 18.4b . where tori inside the 
island can be analytically continued in the complex plane to a manifold escaping to infinity. 

^strictly speaking, outgoing (Sicgcrt) boundary conditions 1471 need to be employed in Eq. 1 18.491 1 in order to properly 
describe the decay process from the well. Those outgoing boundary conditions involve, in addition, an exponential increase of 
the wavefunction's amplitude with increasing distance from the well, which is not taken into account in Eq. 118.4911 assuming 
that the tunncUng rate from the well is comparatively weak. 
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Figure 8.4: Modeling of the direct coupling between the edge of the regular region to the 
chaotic one by a potential barrier separating a potential well from an "open" region, (a) 
Sketch of the potential, (b) Corresponding phase space portrait. 



In that case, Eq. fl8.5ip can be applied provided the tunneling action St is taken as the 
imaginary part of the action integral on a path joining the interior to the exterior of the 
island on this analytical continuation. As a last approximation, one may assume that the 
transiti on t o the "open" part is extremely sharp once the separatrix is crossed. In the model 
of Fig. 18.4b . this amounts to assume a very rapid decrease of the potential, in which case 
one may replace by qm in the tunneling action Eq. fl8.50p . In an actual calculation of the 
direct tunneling rate for a regular state at the edge of the chaos boundary (which we can 
reliably describe only coming from the interior of the island) this amounts to consider in the 
same way that the action S{In, Ic) [E<1- fl8.48p ] provides a good approximation to St- Under 
this hypothesis, one obtains for the coupling to the chaotic Hilbert space the prediction 

., = ^^exp^-2^^J (8.52) 

(see IH] where this computation was proposed with the slightly different language of complex 
time trajectories). 

This "potential-barrier" picture of the direct tunneling is in essence what is behind the 
approach of Podolskiy and Narimanov [19] (though their treatment of the problem is a bit 
mor e soph isticated). Its main virtue is that, beyond the quantized action In = h{n + 1/2), 
Eq. f l8.52p relies only on simple characteristics of the integrable regions : its a rea 2? ?/^ and the 
frequency Qn of the torus In- One needs to keep in mind, however, that Eq. f l8.52p implicitly 
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assum es th at the direct tunnehng mechanism corresponds in some sense to the phase portrait 
of Fig. [8l4b . which in general cannot be justified within any controlled approximation scheme. 



8.3.4 Regular-to-chaotic coupling via a nonlinear resonance 

Another possible, and presumably more realistic, approach to e valua te the direct coupling 
parameter Vn can be obtained assuming that the effective model (18.151) describes the vicinity 
of a resonance not only inside the regular region, but also within the chaotic sea in the 
near vicinity of the island. In this perspective, the model one has in mind for the chaotic 
region follows the spirit of Chirikov's overlapping criterion for the transition to chaos [491 [39] : 
resonances still provide the couplings between quasi-modes, but in the chaotic region these 
couplings become strong enough to completely mix the states. Near the regular-chaos edge, 
the transition between modes inside and outside the regular region is still dominated by one 
or several r:s resonances which might be within or possibly outside the regular island. 

As an illustration, let us consider the simple case where a single nonlinear r:s resonance is 
responsible for all couplings, both within the island and fro m the island's edge to the chaotic 
region. Keeping in mind the discussion in Section 18.2.21 and in Ref . [28] , we assume here 
that the couplings induced by the r:s resonance are dominantly described by the lowest non- 
vanishing Fourier component Vi of the perturbation, i.e. by the matrix elements Vr!^~^^^ = 
{n + r\Hcs\n), and set the phase 0i to zero without loss of generality. 

The structure of the effective Hamiltonian that describes the coupling of the ground state 
Eq to the chaotic sea is, in that case, given by 



( 



Eo Vr 
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y r:s 
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T/[(fcc+l)r 




(8.53) 



where Ekr = E^r — ^Ir-.slkr are the eigenenergies of the unperturbed Hamiltonian T-Lq in the 
co-rotating frame, and the chaotic part (the square in the lower right corner) consists of a 
full sub-block with equally strong couplings between all basis states with actions beyond the 
outermost invariant torus of the islands. In this example the last state within the island 

connected to the ground state is the quasi-mode n = kcV and vl = (W:^^+'^''l)ViVch, with 
Nch the number of states with a given par ity in the chaotic Hilbert space. In a more general 
situation, many resonances may, as in Eq. (18.371) . be involved in connecting the ground stated 
to some |n) at the edge of the island, but we would still write the variance of the matrix 

elements providing the last coupling to the chaotic region as the ratio (W:^+'^'^)ViVch for 
some r:s resonance near the regular-chaos edge. 
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We stress here that it is necessary in this approach to have an exphcit access to the number 
of states A^ch in the chaotic Hilbert space. This is obtained quite trivially for quantum maps, 
such as the kicked rotor we shall consider in section [831 when h = 2tt/N with integer A^. 
In that case is the total number of states in the full Hilbert space and Nch/N represents 
the relative area of the chaotic region in phase space. For a two dimensional conservative 
Hamiltonian, on the other hand, N and thus N^^ are related to the Thouless energy (see 
e.g. the section 2.1.4 of |50]), but this provides only an energy scale rather than a precise 
number, and a more detailed discussion is required to be quantitative. For the sake of clarity, 
we shall in the following limit our discussion to the simpler case that is known. 



8.3.5 Theory of chaos-assisted tunneling 



Let us consider now the effect of the chaotic block on the tunneling process. Eliminating 
intermediate states within the regular island leads for the effective Hamiltonian a matrix of 
the form 



H. 



eff 



V 














TT± 



(8.54) 



for each symmetry class. In the simplest case Eq. f l8.53p where a single r:s resonance needs 
to be considered, the effective coupling matrix element between the ground state and the 
chaos block (H^) is given by 



k=l 



yr-.s 



Eq — E^r + kshw 



.55) 



where En are the unperturbed energies fl8.2Qp of Hcs and \kcr) represents the highest unper- 
turbed state that is connected by the r:s resonance to the ground state and located within 
the island (i.e., Ik^r < h < I(kc+i)r)- More generally, Ves can be expressed in terms of the 
couplings associated with the various resonances that contribute to the t ransitions within 
the island and at the regular-chaos edge. The form of this expression (18.551) already provides 
an explanation for the appearance of plateau-like structures in the tunneling rates. Indeed, 
decreasing h leads to discontinous increments of the maximal number of couplings through 
Eq. (I8.40p and hence to step-like reductions of the effective matrix element Kfj, while in be- 
tween such steps V^ff varies smoothly through the action dependence of the coupling matrix 

(kr) 

elements Vr-.s , provided accidential near-degeneracies in the energy denominators do not 
occur. 

In the simplest possible approximation, which follows the lines of Refs. [121 [17], we neglect 
the effect of partial barriers in the chaotic part of the phase space [H] and assume that the 
chaos block (H^) is adequately modeled by a random Hermitian matrix from the Gaussian 
orthogonal ensemble (GOE). After a pre-diagonalization of (H^), yielding the eigenstates 
(p^ and eigenenergies we can perturbatively express the shifts of the symmetric and 
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antisymmetric ground state energies by 



= Eo + Yl ' (8-56) 



with Ugg_|_ = Vefr(/cr|0^) Performing the random matrix average for the eigenvectors, we 
obtain that {{\{kr\(pf)\'^)) ~ 1/Nch for all j = 1 . . . Nch, which simply expresses the fact that 
none of the basis states is distinguished within the chaotic block (H^). As a consequence, 
the variance of the t'eff±'s is independent of j and equal to v'^^ = V^g/Nch- 

As was shown in Ref. [17], the random matrix average over the eigenvalues £^ gives rise 
to a Cauchy distribution for the shifts of the ground state energies, and consequently also 
for the splittings 

AEo = \E+-E^\ (8.57) 

between the symmetric and the antisymmetric ground state energy. For the latter, we 
specifically obtain the probability distribution 



P{AEo) = (8.58) 



with 



AEo = (8.59) 

^ch 

where A^h denotes the mean level spacing in the chaos at energy Eq. This distribution is, 
strictly speaking, valid only for AEq <^ v^s and exhibits a cutoff at AEq ~ 2vcs, which 
ensures that the statistical expectation value {AEq) = xP{x)dx does not diverge. 

Since tunneling rates and their parametric variations are typically studied on a logarithmic 
scale [i.e., log(A -Eo) rather than AEq is plotted vs. 1/h], the relevant quantity to be calculated 
from Eq. fl8.58p and compared to quantum data is not the mean value (AEq), but rather the 
average of the logarithm of AE^. We therefore define our "average" level splitting {AEo)g 
as the geometric mean of AEq, i.e. 

(A^o), = exp[(ln(AEo))] (8.60) 

and obtain as result the scale defined in Eq. fl8.59p . 



{AEo)g = AEo . (8.61) 



This expression further simplifies for our specific case of periodically driven systems, 
where the time evolution operator U is modeled by the dynamics under the effective Hamil- 
tonian f l8.54p over one period r. In this case, the chaotic eigenphases E^r/h are uniformly 
distributed in the interval [0,27r[. We therefore obtain 



(8.62) 
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for the mean level spacing near Eq. This yields 



{A^o),^-^{AEo),= ^-^j (8.63) 

for the geometric mean of the ground state's eigenphase splitting. Note that this final result 
does not depend on the number A'ch of chaotic states within the sub-block (H^)', as long as 
this number is sufficiently large to justify the validity of the Cauchy distribution f l8.58p (see 
Ref. [IZ]), the geometric mean of the eigenphase splitting is essentially given by the square 
of the effective coupling Ves from the ground state to the chaos. 

The distribution (18.58^ also permits the calculation of the logarithmic variance of the 
eigenphase splitting: we obtain 



([ln(A^o)-(ln(Av.o))]'> = y. (8.64) 

This universal result predicts that the actual splittings may be enhanced or reduced compared 
to (A(y9o)g by factors of the order of exp(7r/2) ~ 4.8, independently of the values of h and 
external parameters. Indeed, as was discussed in Ref. [32], short-range fluctuations of the 
splittings, arising at small variations of h, are well characterized by the standard deviation 
that is associated with Eq. f l8.64p . 

It is interesting to note that the expressio n (I8.63P for the (geometric) mean level spacing is 
quantitatively identical with the expression f l8.45p for the mean escape rate from the regular 
island to the chaotic sea derived in Ref. [22] using Fermi's golden rule. This seems surprising 
as two different nonclassical processes, namely Rabi oscillations between equivalent islands 
and the decay from an island within an open system, underly these expressions. In one- 
dimensional single-barrier tunneling problems, these two processes would indeed give rise to 
substantially different rates; in Eq. (18. ip . to be more precise, the imaginary action integral in 
the exponent would have to be multiplied by two in order to obtain the corresponding expres- 
sion for the decay rate (and the overall prefactor in front of the exponential function should 
be divided by two, which is not important here). The situation is a bit different, however, 
in our case of dynamical tunneling in mixed regular-chaotic systems. In such systems, level 
splittings between two equivalent regular islands involve two identical dynamical tunneling 
processes between the islands and the chaotic sea (namely one process for each island), while 
the decay into the chaotic sea involves only one such process, with, however, the square of 
the corresponding (exponentially suppressed) coupling coef ficien t. This explains from our 
point of view the equivalence of the expressions fl8.63p and (18.450 . 

We finally remark that the generalization of the expression for the mean splittings to multi- 
resona nce pr ocesses is straightforward and amounts to replacing the product of admixtures 
in E q. (I8.55P by a product involving several resonances subsequently, in close analogy with 
Eq. (I8.4ip . In fact, the multi-resonance expression fl8.4ip can be directly used in this context 

replacing the "direct" splittings Ay^i"^ by (K/!s/^''''^V//i)^ where the rf-.Sf resonance is the 
one that induces the final coupling step to the chaotic sea (provided In < Ic < In+rj holds 

for the corresponding action variables; otherwise we would set Aipn^ = 0). This expression 
represents the basic formula that is used in the semiclassical calculations of the splittings in 
the kicked rotor model, to be discussed below. 
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8.3.6 The role of partial barriers in the chaotic domain 

In the previous section, we assumed a perfectly homogeneous structure of the Hamiltonian 
outside the outermost invariant torus, which allowed us to make a simple random-matrix 
ansatz for the chaotic block. This assumption hardly ever corresponds to reality. As was 
shown in Refs. [511 [TTl [12] for the quartic oscillator, the chaotic part of the phase space 
is, in general, divided into several subregions which are weakly connected to each other 
through partial transport barriers for the classical flux (see, e.g.. Fig. 8 in Ref. [I2]) This 
substructure of the chaotic phase space (which is generally not visible in a Poincare surface 
of section) is particularly pronounced in the immediate vicinity of a regular island, where 
a dense hierarchical sequence of partial barriers formed by broken invariant tori and island 
chains is accumulating [521 [531 El] • 

In the corresponding quantum system, such partial barriers may play the role of "true" 
tunneling barriers in the same spirit as invariant classical tori. This will be the case if the 
phase space area AW that is exchanged across such a partial barrier within one classical 
iteration is much smaller than Planck's constant 2'irh [55], while in the opposite limit AW ^ 
27rh the classical partial barrier appears completely transparent in the quantum systen:0. 
Consequently, the "sticky" hierarchical region around a regular island acts, for not extremely 
small values of h, as a dynamical tunneling area and thereby extends the effective "quantum" 
size of the island in phase space. As a matter of fact, this leads to the formation of localized 
states (also called "beach" states in the literature [T3]) which are supported by this sticky 
phase space region in the surrounding of the regular island [56] (see Fig. 18.6b ). 

An immediate consequence of the presence of such partial barriers for resonance-assisted 
tunneling is the fact that the critical action variab le Ic deflning the number kc of resonance- 
assisted steps within the island according to Eq. ( 18.55^ should not be determined from the 
outermost invariant torus of the island, but rather from the outermost partial barrier that 
acts like an invariant torus in the quantum system. We flnd that this outermost quantum 
barrier is, for not extremely small values of h, generally formed by the stable and unstable 
manifolds that emerge from the hyperbolic periodic points associated with a low-order non- 
linear r:s resonance. These manifolds are constructed until their flrst intersection points in 
between two adjacent periodic points, and iterated r — 1 times (or r/2 — 1 times in the case 
of period-doubling of the island chain due to discrete symmetries), such as to form a closed 
artiflcial boundary around the island in phase spac43 As shown in Fig. 18. 5[ one further 
iteration maps this boundary onto itself, except for a small piece that develops a loop-like 
deformation. The phase space area enclosed between the original and the iterated boundary 
precisely deflnes the classical flux AW exchanged across this boundary within one iteration 
of the map [52l [53] . 

The example in Fig. 18.51 shows a boundary that arises from the inner stable and unstable 
manifolds (i.e. the ones that would, in a near-integrable system, form the inner separatrix 
structure) emerging from the unstable periodic points of a 4:1 resonance (which otherwise 
is not visible in the Poincare section) in the kicked rotor system. Judging from the size of 
the flux area AW, this boundary should represent the relevant quantum chaos border for 
the tunneling processes that are discussed in the following section. We clearly see that it 



^More precisely, the authors of Ref. 1551 claim that AW has to be compared with nh in order to find out whether or not a 
given partial barrier is transparent in the quantum system. 

^''This construction is also made in order to obtain the phase space areas S^,s that are enclosed by the outer and inner 
separatrix structures of an r:s resonance, and that are needed in order to compute the mean action variable I^.s and the 
coupling strength Vr-.s of the resonance according to Eqs. 118.261 1 and US. 271 1. 
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encloses a non-negligible part of the chaotic classical phase space, which includes a prominent 
10:3 resonance that, consequently, needs to be taken into account for the coupling process 
between the regular island and the chaotic sea. Thereby, we naturally arrive at multi-step 
coupling processes across a sequence of several resonances, which would have to be computed 
for a reliable prediction of the tunneling rates in the semiclassical regime. 




e 



Figure 8.5: Classical phase space of the kicked rotor at = 3.4 in approximate action- angle 
variables (1,0). The thick solid line shows the location of the effective quantum boundary 
of the central island for the values of Planck's constant that are considered in Section 18.41 
This effective boundary is constructed from segments of the stable and unstable manifolds 
that emerge from the hyperbolic periodic points of the 4:1 resonance at ^ ~ 7r/6 and 57r/6, 
respectively. Those segments were computed until the symmetry axis at 6 = 7r/2 and then 
iterated three times under the kicked rotor map, yielding the thick solid line. A further 
iteration of this boundary maps it onto itself, except for the piece between 6 ^ O.IStt and 
6 = O.Svr which is replaced by the lighter curve. The phase space area that is enclosed 
between the original (dark) and the iterated (light) boundary defines the classical flux that 
is exchanged across this boundary within one iteration of the map. The dashed line shows, 
in comparison, the actual classical chaos border defined by the outermost invariant torus of 
the island. 



8.4 Application to the kicked rotor 



8.4.1 Tunneling in the kicked rotor 

To demonstrate the validity of our approach, we apply it to the "kicked rotor" model, which 
is described by the Hamiltonian 

oo 

H{p,q,t)=p'^/2-K 6{t -n) cos q. (8.65) 



n=— oo 
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The classical dynamics of this system is described by the "standard map" {p, q) i— )■ {p', q') 
with 

p' = p — K sin q (8.66) 

q' = q + p, (8.67) 

which generates the stroboscopic Poincare section at times immediately before the kick. The 
phase space of the kicked rotor is 27r periodic in position q and momentum p, and exhibits, 
for not too large perturbation strengths K < 4, a region of bounded regular motion centered 
around {p,q) = (0,0). 

The quantum dynamics of the kicked rotor is described by the associated time evolution 
operator 

[/ = exp I — — — ) exp ( — -^i^cosg ) (8.68) 
\ h2 J \ h J 

which contains two unitary operators that describe the effect of the kick and the propagation 
in between two kicks, respectively {p and q denote the position and momentum operators). 

Because of the classical periodicity in both p and g, we can consider tunneling between 
the main regular island centered around (p, q) = (0, 0) and its counterparts that are shifted 
by integer multiples of 27i along the momentum axis or along the g-axis. To mimic a double- 
well configuration, we will restrict the boundary conditions for the eigenstates of U and 
consider tunneling between two islands centered around (i) {p, q) = (0, 0) and (27r, 0) or 
around (ii) {p,q) = (0,0) and (0,27r). The effective parity that allows to discriminate the 

eigenphases ip^ of U manifests as 'ipniP + 2'/r) = ±ip^{p) for the corresponding eigenstates 
in momentum representation in case (i) and as ipnil + 2vr) = ±ip^{q) for the eigenstates 
in position representation in case (ii). In both cases, tunneling will be characterized by the 
splitting 

A<^„ = |<^+-<^-| . (8.69) 

Numerically, it can be convenient to deal a finite Hilbert space of (even) size A^, and this 
can be obtained provided the two phase space translation operators Ti = exp{27iip/h) and 
T2 = exp{—2'7iiq/h) commute, which is the case if we choose h = 271 /N [571 158] . 



8.4.2 Numerical computation of the eigenphase splittings 



Figure 18.61 shows the eigenphase splittings Ac^jq [see Eq. fl8.69p ] in case (ii) (tunneling in 
position) for the local "ground state" (n = 0) in the central island of the kicked rotor, i.e. for 
the state that is most strongly localized around the center of the island, a.t K = 2.28. While 
on average these splittings decrease exponentially with 1/h, significant fluctuations arise on 
top of t hat exponential decrease. In particular, large spikes are visible. As illustrated in 
Fig. 18.6b . they can be related to the crossing of "excited states" within the island, which are 
coupled to the ground state by a classical resonance. Fig. 18.6b shows the Husimi distribution 
of the relevant states involved, demonstrating that the coupling process is most effective when 
the states are symmetrically located on each side of the classical resonance. To illustrate 
that the influence of the resonances inside the regular islands is actually independent of the 
details of the chaotic regime (a major feature of th e resonance-assisted and chaos-assisted 
tunneling schemes) we furthermore plot in Fig. 18.71 a comparison between the splittings for 
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Figure 8.6: The upper panel shows the quasi-energy sphttings in AEq = h(po in the kicked 
rotor model for K = 2.28, corresponding to a 4 : 1 classical resonance. Here we are concerned 
with a tunneling in q. The middle panel shows the quasi-energy spectrum E^ = hip^ of U 
where only states with a significant overlap with a coherent state localized around (p, q) = 
(0, 0) have been retained. The horizontal arrow a) marks the central state doublet (n = 0, 
not resolved at that scale) and the arrow b) indicates the third excited state localized in 
the island (n = 4). Their Husimi distribution superimposed with the Poincare surface of 
section are shown in the lower panel in order to illustrate the clear correspondence between 
the classical and the quantum resonance. As 1/h increases, the crossing of the doublet by 
the resonant state provokes the large and wide spike indicated by the vertical arrow in the 
upper panel. 
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the cases (i) and (ii) of tunneling in p and q direction, respectively. The (rough) matching 
between the dominant spikes of fluctuations in both cases confirms that tunneling outside 
one island is mainly isotropic in phase space. 
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Figure 8.7: Comparison between the splitting Aipo for the kicked rotor at i^' = 2 for (i) the 
case of tunneling between two islands centered at (p, q) = (0, 0) and at (p, q) = (27r, O), and 
(ii) the case of tunneling between two islands centered at (p, q) = (0, 0) and at {p, q) = (0, 27i). 
In this latter case, we restricted ourselves to even integer values of = 27T/h for which 
the splittings can be computed by diagonalizing finite N x N matrices beyond the double 
precision. 

From now on, we will consider case (i) only (tunneling in momentum) with = 27r/h 
being a n int eger , corr esponding to the number of Planck cells that fit into one Bloch cell. 
Figures 18.91 and 18.101 show the eigenphase splittings Aipo of th e isla nd's ground state for 
K = 2.6, 2.8, 3.0 (Fig. EB as well as for K = 3.2, 3.4, 3.6 (Fig. [8A0|) . As in Refs. [SH E^, 
these sphttings were calculated with a diagonalization routine for complex matrices that is 
based on the GMP multiple precision library jSH], in order to obtain accurate eigenvalue 
differences below the ordinary machine precision limit. 



8.4.3 Semiclassical calculations 

The role played by the nonlinear resonances in the tunneling mechanism is made explicit by 
comparing these numerically calculated splittings with semiclassical predictions based on the 
most relevant resonances that are encountered in phase space. In practice, we took those r:s 
resonances into account that exhibit the smallest possible values of r and s for the winding 
numbers s/r under consideration. In all of the considered cases, th e "qu antum boundary" of 
the regular island, which determines the value of kc through Eq. fl8.40p . was defined by the 
partial barrier that results from the intersections of the inner stable and unstab le m anifolds 
associated with the hyperbolic periodic points of the 4:1 resonance (see also Fig. 18. 5p . While 
this partial barrier lies rather close to the classical chaos border of the isla nd fo r K = 2.6 
(Fig. 18. 9j) . it encloses an appreciable part of the chaotic sea for K = 3.6 (Fig. I8.10p including 
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some relevant nonlinear resonances. 



Pure resonance-assisted tunneling 

We stress that the semiclassical calculations shown in Figs. 18.91 and 18.101 involve a few dif- 
ferences as compared to some of our previous publications [SL .32^ 33., 34J . To start with, (I) 
the action dependence of the coupling coefficients associated with the resonances has been 
included [see Eq. fl8.36p ]. Furthermore, the unperturbed energy differences E n — E n+kr of the 
quasi-modes are not computed via the quadratic pendulum approximation f l8.24p . Instead, 
as illustrated on Fig. 18. 8[ (II) a global parabolic fit to the action dependence of the frequency 
Q = fi(J) was applied on the basis of the classically computed values of the resonant actions 
Ij.;s and their frequencies Q{Ir:s) = {s/r)u, for a sequence of resonances with not too large 
r and s. Theses two modifications significantly improve the reproduction of individual peak 
structures in the tunneling rates. 

unperturbed energy H (I) 




action variable I action variable I 



Figure 8.8: Unperturbed energies and oscillation frequencies within the regular island of 
the kicked rotor at K = 3.5. The left panel shows the action dependence of the oscillation 
frequencies as computed from a quadratic fit to individual nonlinear resonances (crosses), 
which results in the expression n{I) = 2.41886 - 0.790561J + 0.235191/^ (sohd line). The 
upper right panel compares the unperturbed energies resulting from the integration of this 
quadratic expression (solid line) with the unperturbed energies used in Ref. [35] that were 
obtained by analyzing a dense set of quasi-periodic trajectories within the regular island 
(dashed line) [60]. This latter approach yields n{I) = 2.41740 - 0.952917/ + 1.00151/^ - 
1.00153/^ + 0.368829/^, which essentially constitutes the definition of the fictitious integrable 
system used in Ref. [35]. Although the difference between these two approaches is rather 
small as shown in the lower right pa nel, it plays a significant role for the tunneling rates in 
the deep semiclassical limit (see Fig. I8.12p . 

With these improvements (I) and (II), we generally find that the quantum splittings are 
quite well reproduced by our simple semiclassical theory based on nonlinear resonances. 
In particular, the location and height of prominent plateau structures and peaks in the 
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tunneling rates can, in almost all cases, be quantitatively reproduced through resonance- 
assisted tunneling. The additional fluctuations of the splittings on a small scale of N, 
however, cannot be accounted for by our approach as they arise from the details of the 
eigenspectrum in the chaotic block of the Hamiltonian. Their average size, however, seems in 
good agreement with the universal prediction (I8.64p for the variance of eigenphase splittings 
in chaos-assisted tunneling. 

Note that there is a general tendency of the semiclassical theory to overestimate the 
exact splittings wherever the latter encounter local minima. We attribute those minima 
to the occurrence of destructive interferences between different pathways that conne ct th e 
ground state to a given excited state \n). As pointed out in the discussion of Eq. (I8.4ip . 
such destructive interferences are not yet accounted for in our present implementation of 
resonance-assisted tunneling. Their inclusion woul d req uire to take into account the phases 
4>k associated with individual resonances [see Eq. (18.131) ]. the discussion of which is beyond 
the scope of this article. 



Resonance-assisted tunneling at a bifurcation 



As a last comment, we note that the influence of a nonlinear resonance on tunneling pro- 
cesses in mixed systems may persist even if that resonance is not at all manifested in the 
classical phase space. This is precisely the case at the value of the perturbation parameter 
at which this resonance is bifurcating from the center of the island, i.e., at which the central 
fixpoint of the island exhibits a rational winding number corresponding to the resonance 
under consideration. A prominent example in the kicked rotor model is found a.t K = 2 
where the central fix point has the winding number 0.25. Indeed, = 2 is exactly the 
critical value where two periodic orbits of period 4 coalesce in the center, both separately 
coming from the complex phase space (their action being strictly negative for K < 2) and 
then becoming real and distinct (one stable and one unstable) for K > 2. 

Following the normal-form arguments in Section I8.2.3[ the scaling of the "classical size" 
of the resonance with the perturbation parameter K is, in lowest order, provided by the 
effective pendulum matrix element Vr:s{K) = Vi[Ir:s{K)Y''^^ where Ir-AK) represents the 
dependence of the resonant action on K. However, as can be seen in Eq. (18.361) . the associated 
coupling matrix elements that affect tunneling only depend on the prefactor {ii, which ought 
to be a well-behaved function of K showing no singular behaviour at the bifurcation point. 
Therefore, we have to conclude that these matrix elements remain finite even for K <2. 

This expectation is confirmed in Fig. 18.111 which shows the tunneling rates in the kicked 
rotor for K = 2. For this value, the only major resonance that is manifested in the classical 
phase space is the 10:2 resonance whose island chain is located near the boundary of the 
main island. Taking into account this resonance alone (as well as combining it with other 
resonances of higher order, the result of which is not shown in Fig. I8.1ip apparently leads to 
a very strong underestimation of the quantum eigenphase sphttings (in contrast to Ref. [31] 
where the action dependence of the coupling matrix elements was not properly incorporated). 
But once we take into account the 4:1 resonance and compute its corresponding classical 
parameters from the classical phase space ai K = 2.001, we obtain a good reproduction of 
the quantum splittings. 
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Figure 8.9: Quantum and semiclassical splittings in the kicked rotor model for K = 2.6 (left 
column) K = 2.8 (central column), and K = 3 (right column). The upper and middle panels 
show the classical phase space in the original phase space variables p and q, with the thick 
curve marking the effective quantum boundary of the island, and in approximate action- 
angle variables / and 6. The lower panels display the quantum and semiclassical eigenphase 
splittings (dots and solid lines, respectively) of the ground state in the central regular island. 
For the semiclassical splittings, we used the 14:4 and 18:5 resonances for K = 2.6, the 10:3 
and 14:4 resonances for K = 2.8, and the 10:3, 14:4, 16:5, and 22:7 resonances for K = 3. 
As pointed out in the text, the splittings were computed with a generalization of the multi- 
resonance expression (18.411) to mixed systems, using (I) the corrected action dependence 
(18.361) of the matrix elements and (II) unperturbed energies that were determined from a 
global parabolic fit of fi(/), and (III) computing the coupling to t he cha otic domain via the 
outermost nonlinear resonance, as explained at the end of section I8.3.5[ The discontinuous 
steps in the semiclassi cal sp littings are actually induced by the discontinuity of the integer 
part in the expression (I8.40p for the maximal number kc of couplings within the island. 
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Figure 8.10: Same as Fig. 18.91 for K = 3.2, K = 3.4, and K = 3.6. For the semiclassical 
splittings, we used the 6:2 and 10:3 resonances for K = 3.2, the 6:2, 10:3, and 14:5 resonances 
for K = 3.4, and the 6:2 and 8:3 resonances for K = 3.6. 



Combination with direct regular-to-chaotic couplings 



In the semiclassical predictions shown in Figs. 18.91 18.101 and 18.111 we have, following the 
discussion in section 18.3.41 assumed that the resonance-assisted mechanism was providing 
the relevant couplings not only within the regular island, but als o (III) from the edge of 
the regular island to the chaotic sea. As pointed out in section I8.3.2[ the regular-chaos 
boundary is, however, the place where approximation schemes are not controlled any longer. 
It is therefore useful to compare the results obtained in this way with those derived from 
the more precise evaluation of the direct regular-chaotic couplings (Illb) computed with 

Eq. ( ]8.45p through a numerical application of the quantum evolution operator U. 



Figure 18.121 shows the resulting comparison for the eigenphase splittings of the quantum 
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Figure 8.11: Resonance-assisted tunneling in the kicked rotor at a bifurcation. The left 
panel shows the classical phase space for K = 2, which contains a prominent 10:2 resonance 
close to the border of the island, and the right panel displays the corresponding quantum 
eigenphase splittings (dots). Semiclassical calculations of the splittings (solid and dashed 
lines, using the same levels of approximation as in Figs. 18.91 and I8.10p were carried out using 
the 10:2 resonance only (dashed line) as well as a combination of the 4:1 resonance and the 
10:2 resonance, the former emerging at the island's center right at K = 2. The parameter vi 
associated with that 4:1 resonance [see Eq. fl8.36p ] was determined from the classical phase 
space at K = 2.001. 

kicked rotor at K = 3.5 [35] (see Fig. 18.31 for the corresponding classical phase space). 
In addition to the quantum and semiclassical splittings obtained in the same way as in 
Figs. 18.91 and 18.101 (black and red curves, respectively), two additional curves are shown. 
The green one corresponds to a fully semiclassical calculation for which the regular-to-chaotic 
couplings are evaluated by the reson anc e-assi sted mechanism. In contrast to the red curves 
and to the calculations in Figs. 18.91 and [8.10[ however, the action dependence Hq{I) of the 
unperturbed energies within the island was not obtained by a fit to several relevant resonances 
as described above, but rather (lib) by computing the action and the rotation number for 
a dense set of trajectories within the island, from which the energies are deduced via the 
relation fl{I) = dHQ{I)/dI [60] (see, e.g.. Fig. 11 of Ref. [H] and the associated text for a 
detailed discussion). The blue curve implements, in addition to these improved energies, an 
evaluation of the direct regular-chaotic coupling (Illb) using the approach of Backer et al. 
[22] [see Eq. (18.451) ]: it actually corresponds to the curve published in Fig. 3 of Ref. [35] . 

What is observed in Fig. 18. 121 is that, although there is good agreement between the three 
theoretical curves and the numerical one for a large range of iV = 27r/h, some significant 
deviations arise in the range 300 < N < 450. Remarkably, however, the most striking dis- 
crepancies that are encountered for the red curve (i.e., for the full semiclassical calculation 
without the improved energies) are essentially cured once the improved energies are imple- 
mented (green curve). Adding the semi-numerical evaluation of the direct regular-to-chaotic 
coupling further improves the prediction, but to a significantly lesser degree. 



Figure \8M compares Hq{I) computed for the kicked rotor at K = 3.5 by the two methods 
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Figure 8. 12: Quantum and semiclassical splittings in the kic ked rotor model for K = 3.5 
(see Fig. 18.31 for the classical phase space). As in Figs. 18.91 and 18.101 the dots and solid 
lines represent, respectively, the quantum splittings and the semiclassical prediction. In the 
upper left panel, this prediction is based on our approach, where we take into account the 
6:2, 8:3, 10:3, and 14:5 resonances. The upper right panel shows the same calculation, except 
that a more refined evaluation of the unperturbed energies (and their action dependence) 
is used (lib) (see text and Fig. 18. 8p . The lower left panel shows the prediction that is 
obtained with these improved energies and with an evaluation of the direct regular-chaotic 
coupling (Illb) using the fictitious integrable system approach of Backer et al. [22] [which 
amounts to applying the qu antum kicked rotor map onto the corresponding eigenstate of the 
integrable system, see Eq. f l8.45p ]. Finally, the lower right panel shows the prediction that 
would be obtained with the 6:2 resonance according to the approach outlined in Ref. j3l], i.e. 
neglecting the action dependence of the coupling matrix elements, neglecting the occurence 
of partial barriers in the chaotic domain, and making a simple quadratic expansion of the 
unperturbed energies in the vicinity of the resonance. 



under consideration, namely the fitt ing approach using the most relevant resonances (which 
is illustrated in the left panel of Fig. 18. 8p and the more refined approach based on a dense set 
of trajectories [60]. Apparently, both approaches yield nearly identical energies, with relative 
differences being well below one percent (see the lower right panel of Fig. 18. 8p . This underlines 
that small imprecisions in the prediction of the unperturbed eigenenergies of regular quasi- 
modes may, under special conditions, lead to dramatic over- or underestimations of the 
tunneling rates in mixed systems. In the particular case considered here, it seems to be the 
transition across the 8:3 resonance, located rather close to the center of the island, which is 
not properly described using the more approximate energies. This gives rise to a horizontal 
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shift of the predicted sphttings towards smaller A^, as is clearly seen in Fig. 18.121 This issue 
obviously requires further investigations. It does, however, not seem to put into question the 
principal conclusion that the resonance-assisted coupling mechanism provides an accurate 
approach to evaluate regular-to-chaotic tunneling rates in a purely semiclassical manner. 



8.5 Conclusion 



In summary, we have provided a comprehensive description of the theory of resonance- 
assisted tunneling in mixed regular-chaotic systems. This description is partly based on 
previous publications of ours [27l EH [HI |32l |33], but contains also some new aspects and 
significant improvements especially concerning the determination of the matrix elements as- 
sociated with resonance-induced couplings. Moreover, partial barriers in the chaotic domain 
are now incorporated into the scheme of resonance-assisted tunneling, which generally gives 
rises to an effective enhancement of the size of the regular islands under consideration. In 
practice, we find that the most relevant partial barriers for tunneling are constituted by 
combinations of stable and unstable manifolds that are associated with hyperbolic periodic 
points of low-order nonlinear resonances within the chaotic domain. The application of this 
approach to the kicked rotor model yields rather good agreement between the "exact" eigen- 
phase splittings of states that are localized in the center of the main regular island, and their 
semiclassical predictions based on nonlinear resonances. 

The main message that we intend to communicate here is that resonance-assisted tun- 
neling not only allows one to understand the origin of plateaus and peak structures, in the 
tunneling rates. It also provides a simple, readily implementable scheme to quantitatively 
predict the appearance of such structures on the basis of purely classical information. In 
practice, one needs for this purpose to identify the relevant resonances in the regular phase 
space region under consideration, to find their stable and unstable fixpoints in the Poincare 
surface of section, to compute their stability indices and the areas enclosed by their separa- 
trix structures, respectively by their stable and unstable manifolds, and finally to compute 
the flux enclosed by the turnstiles in order to determine effective "quantum chaos border" 
of the island. Even for a simple model like the kicked rotor, this programme requires much 
less numerical effort than a quantum calculation of the tunneling rates. As we demonstrated 
for the kicked rotor, it provides, on the other hand, a reproduction of the quantum tunnel- 
ing rates which is extremely satisfactory from a quantitative point of view. This numerical 
accuracy requires, however, a careful evaluation of the various classical parameters entering 
into the semiclassical calculation of the splittings. This includes the action dependence of 
the resonance-Hamiltonian couplings, the effective size of the regular island, and the evalu- 
ation of the unperturbed energies within the island. It turns out in particular that in some 
circumstances, e.g. a.t K = 3.5 in the range 300 < N < 450, very small imprecision in the 
determination of the unperturbed eigenenergies may significantly affect the accuracy of the 
semiclassical predictions. 

We expect that the framework of resonance-assisted tunneling can be generalized to sys- 
tems with more than two effective degrees of freedom, although the identification of im- 
portant resonances might become more involved in such systems and resonance-assisted 
(quantum) tunneling might compete there with (classical) Arnol'd diffusion in the deep 
semiclassical regime [30]. Another open problem which needs to be addressed in more de- 
tail concerns the role of nonlinear resonances in trajectory-based semiclassical approaches 
to tunneling, put forward by Shudo, Ikeda, and coworkers [361 EZl [38] (see also the corre- 
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sponding chapters in this book), which involve the complexified classical phase space and 
are more rigorous from a formal semiclassical point of view than our approach. We strongly 
believe that nonlinear resonances should leave their characteristic traces in the self-similar 
complex phase space structures that govern tunneling in this framework [5^ [57] . Such an 
insight should significantly contribute to rendering those approaches practicable for more 
complicated systems as well — and underline the semiclassical nature of resonance-assisted 
tunneling. 
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